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Abstract 

We calculate the inclusive production of charged hadrons in pp collisions to next-to-leading order 
(NLO) in the QCD improved parton model using a new set of NLO fragmentation functions for 
charged pions and kaons. We predict transverse-momentum distributions and compare them with 
experimental data from the CERN SppS Collider and the Fermilab Tevatron. 
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1 Introduction 

The inclusive production of single hadrons in hadron- 
hadron, photon-hadron and deep inelastic lepton- 
hadron collisions is an important area to test the 
QCD improved parton model. The inclusive cross sec- 
tion is expressed as a convolution of the parton dis- 
tribution functions, the partonic cross sections and 
the fragmentation functions of quarks and gluons into 
charged or neutral particles. The factorization the- 
orem ensures that parton distribution and fragmenta- 
tion functions are universal functions and that only the 
hard scattering partonic cross sections change when 
different processes are considered. This theory pro- 
vides a rather consistent description of many large- 
momentum-transfer-processes [1] and it is suitable for 
describing the inclusive production of single hadrons 
at large transverse momentum (px) in various reac- 
tions. In this work we study the inclusive production 
of charged pions and kaons in high energy pp collisions. 

Experimental data for charged single particle pro- 
duction come from the CERN ISR pp Collider [2] ; from 
the UA1 [3] and UA2 [4] Collaborations at the SppS 
Collider at CERN; and from the CDF Collaboration 
at the Tevatron [5] . Recently, high statistics data from 
the UA1 MIMI Collaboration [6] have become available 
which has extended the px range to much larger values 
as compared to earlier UA1 analysis. 

In [7], together with B.A.Kniehl, we already pre- 
sented results for inclusive single-charged-hadron and 
singlc-7T° cross sections at the next-to-leading order 
(NLO) [7]. We compared them with experimental data 
from the UA2 and CDF Collaborations. The overall 
agreement concerning the px dependence and the ab- 
solute normalization of the cross sections between the- 
oretical and experimental results was satisfactory even 
for the smaller px region. An important drawback of 
this work, however, was due to the use of leading order 
(LO) parametrizations of fragmentation functions for 
charged pions and kaons [8] . These old parametriza- 
tions had been constructed more than ten years ago 
from fits to the data then available, obtained in low- 
energy e + e~ annihilation experiments and deep inelas- 
tic muon-nucleon scattering. 

Recently a NLO set of parametrizations for frag- 
mentation into charged pions and kaons was ob- 
tained [9] . These parametrizations are generated 
through fits to e + e _ annihilation data taken at y/s = 
29GeV by the TPC Collaboration. They produce 
rather satisfactory fits also for other e + e~ data for 
charged particle production obtained at lower energy 
at DORIS and at higher energy at PETRA, PEP and 
LEP. 



By using these new parametrizations, we are now 
in position to perform a full NLO calculation includ- 
ing NLO parton distribution functions, NLO parton- 
parton hard scattering cross sections and NLO frag- 
mentation functions thereby, removing a serious limi- 
tation of our earlier work [7] . 

Other authors have almost simultaneously pre- 
sented NLO parametrizations of fragmentation func- 
tions for light mesons, i.e. neutral pions [10], eta 
mesons [11] and charged pions [12]. In general these 
parametrizations are obtained through fits to data pro- 
duced by the HERWIG Monte Carlo [13] at fixed Q = 
30 GcV. Parametrizations at lower scales are then ob- 
tained via NLO evolution. The authors of [10] and [11] 
compared their 7r°-fragmcntation functions with data 
from pp collisions in fixed target experiments and the 
ISR, and from pp collisions at the CERN SppS Col- 
lider. In [10] also, several sets of fragmentation func- 
tions for 7T° and rj were obtained at a low scale, com- 
parable to the one used in [8] through simultaneous fits 
to e + e~ annihilation, fixed target and ISR collisions 
and SppS Collider data. 

In this work we follow a different approach. We as- 
sume that the NLO fragmentation function for charged 
pions and kaons are sufficiently constrained by the fit 
to the TPC e + e" annihilation data and the fit to the 
gluon fragmentation function from the OPAL measure- 
ments [14]. Our aim is to verify whether the NLO 
parametrization of [9] gives also a satisfactory account 
of existing pp collider data and whether it improves the 
estimates obtained in our earlier work [7] . As in [7] , we 
assume that the charged particle yield can be very well 
described by the sum of the pion and kaon yield and 
that the inclusive production of other charged particles, 
like p and p and heavier baryons is negligible compared 
to the production of the two lighter mesons. 

This work is organized as follows. In Sect. 2 we 
give a short introduction of the formalism and fix our 
inputs. The results of our calculation and a comparison 
with data from pp colliders are presented in Sect. 3. 
This sections ends with a discussion of the results and 
concluding remarks. 



2 Formalism and Input 

Before presenting the results, we give a brief introduc- 
tion to the NLO formalism. The NLO inclusive cross 
section for the production of a single hadron h in the 
reaction 

P(pi) + P(P2) -> hfa) + X , (1) 
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Fig. 1. Inclusive cross section for production of charged- 
hadrons (h = (h + + h~)/2) as a function of pt for y/s — 
200 GeV and rapidity \y\ < 2.5. The short-dashed, solid 
and long-dashed lines correspond to the full NLO prediction 
for scales fi, M and Mf set equal to Pt/2, pt and 2pT- For 
comparison, the UA1 data [3] taken at the same energy and 
in the same rapidity range are also shown 

is written as: 

a,b,c J 

x F*( Xl ,M 2 ) F^(x2,M 2 ) D h c (x 3 ,Mj) 



irS 



dv 



6(1 -w) 



+ 



2ir 



K ab ^ c (s,v,w;^ 2 ,M 2 ,M 2 f ) 



(2) 



The partonic variables v and w are related to the 
usual s,t,u (s = (p a +p b ) 2 , t = (p a - p h ) 2 and u = 
(pt - Pc) 2 ) as: v = 1 + t/s, w = —u/(s + t). They 
are related to the hadronic variables S = (p\ +P2) 2 , 
T =(pi- p 3 ) 2 and U = (p 2 - p 3 ) 2 by: 



s = xix 2 S , 



t=^T, 

X3 



X3 



(3) 
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The indices a,b,c run over gluons and Np flavours of 
quarks. As in our earlier work [7], we assume Np = 4 
and neglect the influence of the charm-quark threshold. 
FP(xi,M 2 ) and D^(x 3 ,M 2 ) are the usual structure 
and fragmentation functions for partons of type a and c 
respectively inside the proton and the hadron h. They 
depend on x\, x 3 , the fractions of proton momentum 
carried by parton a and the fraction of the momentum 
of parton c carried by hadron h and on the factorization 
scales M and Mf. The additional mass parameter p, 
is the renormalization scale for the strong coupling. 
Finally, d 3 a® b ^ c is the LO partonic cross section for 
the process a + b ^ c + X in O (a 2 (p 2 )). The K ab ^ c 
functions are the NLO corrections to a + b^c + X 
and are taken from the work by Aversa et al. [15]. 

The form of the coefficients K ab ^ c is not unique. 
They depend on the choice of finite corrections fcj (x) 
and dij (i,j = q,g) to structure and fragmentation 
functions. These subtraction terms are accompanied 
by appropriate definitions of the FJ- and We 
choose here the structure functions in the MS scheme. 
Specifically, we take the CTQ2M (MS) set of the pro- 
ton distribution function from [16] which gives a good 
fit to recent structure function measurements at small 
a; [17]. The starting scale of this set is M = 2 GeV 
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Fig. 3. Same as in Fig. 1 for V* = 900 GeV 

and therefore it can be used for prediction of the sin- 
gle particle cross section at fairly small values of p^. 
The a s is computed from the two-loop formula with 
Aj^- = 0.231 GeV as in the fit of the CTQ2M (MS) 
structure function. 

As functions £)■?, we employ the parametrizations 
by Binnewies et al. [9]. They are given separately for 
the average of charged pions and kaons, also in the MS 
scheme. They are extracted from the e + e~ data on in- 
clusive production of pions and kaons at Q — 29 GeV 
and found in agreement with the majority of the e + e~ 
data in the energy range between 5 GeV and 91 GeV. 
The Aj^g = 0.19 GeV chosen for these fits is com- 
patible with the A-jg-j in the proton structure func- 
tion. The starting energy Qq, where fragmentation 
functions have simple parametrizations, is chosen to 
be Qo — \/2GeV. The x and Q 2 dependence is then 
obtained by a NLO evolution. 

The gluon fragmentation function into pions and 
kaons is poorely determined in these fits to e + e~ an- 
nihilation data. This is due to the fact that the gluon 
participates in the process only at NLO. It has an ap- 
preciable impact on the cross section only at very small 
values of x, contributing mainly through the Q 2 evo- 
lution, due to its coupling to the singlet combination 
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Fig. 4. Inclusive cross section for production of hadrons 
h = h + + h~ and kaons, K = K + + K~ as a function of pt 
for y/s — 630 GeV and rapidity range \y\ < 3.0. The data 
obtained by UA1 MIMI Collaboration [6] for production of 
hadrons are also shown 

of quarks. Therefore the gluon parameters are corre- 
lated with the sea quark parameters and are not very 
well constrained by e + e~ data. To fix them indepen- 
dently, the gluon fragmentation into charged particles 
is compared in [9] to the three-jet data of the OPAL 
Collaboration at LEP [16]. 

In an earlier attempt [17] when this last compari- 
son was not made, an equally good description of all 
the other e + e~ data used in [9] could be obtained with 
weaker gluon fragmentation functions. To show the 
gluon effect, we shall give for comparison, the inclusive 
charged particle cross section with the weaker gluon 
distribution of [7] for one of the pp energies. 

Further details of the fragmentations and their 
parametrization as a function of Q 2 can be found in [9] . 



3 Numerical results 

We are now in position to present our numerical results. 
We work in the MS scheme with Nf — 4 active quark 
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Fig. 5. The ratio of inclusive cross sections for charged K, 
K = K + + K~ , over charged hadrons h, h = h + + hT , as 
a function of pt for \fs = 630 GeV and \y\ < 3.0 



flavours. We set the three scales /i, M, Mf equal and 
we vary them between px/2 and 2px- Unless otherwise 
specified the charged hadron h is defined as h = (h + + 
h~)/2, where h 



sums over 



7r and K 



In Figs. 1,2 and 3 we show the inclusive charged 
hadron cross section for p+p — > /i+X at y^I = 200, 500 
and 900 GeV. The rapidity is averaged over the interval 
-2.5<y<2.5. The agreement with the UA1 data [4] 
is best with scales equal to px, except at y/s = 200 GeV 
where the data lie somewhat nearer to the prediction 
with scales equal to px/2. In agreement with the exper- 
imental data, the theoretical curves show the expected 
increase of the high px tail between y/s — 200 GeV and 
y/s = 900 GeV. According to [3] , additional systematic 
errors due to luminosity and acceptance corrections are 
small (±15%). Therefore we can conclude that our pre- 
dictions agree well in shape and normalization with the 
data. The data with small pt have the smallest experi- 
mental errors. Unfortunately below px = 2-3 GeV, our 
predictions cease to be valid: the soft production mech- 
anism takes over while the hard production mechanism 
tends to give too fastly growing results. 

The most recent and best experimental data come 
from the UA1 MIMI Collaboration [6] . They have the 
smallest experimental errors and extend earlier analysis 
of UA1 measurements to px values up to 25 GeV. These 
data are relative to production of charged hadrons h, 
with h = h + + h~ . The theoretical predictions ob- 
tained for the same definition of h are plotted in Fig. 
4 for scales pt/2,pt and 2pT, together with the data 
from [6] . The data and the theoretical predictions are 
for \fs = 630 GeV and are averaged over \y\ < 3.0. 
The agreement with the data from the small up to the 
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Fig. 6. Same as in Fig. 1 for ^fs = 630 GeV and rapidity 
range \y\ < 1.0. The theoretical results are compared with 
the data obtained by the CDF [5] 

highest values of px is excellent. The curve with scale 
Pt/2 fits the data best over the whole px range above 
3 GeV. 

In Fig. 4 we give also our prediction for K produc- 
tion, in case such data become available in the future. 
The ratio to the total charge particle K/h is plotted 
in Fig. 5. Above px = 5 GeV this ratio is approxi- 
mately constant in px with a slight increase towards 
larger px- This ratio is approximately 0.15 and reflects 
the fact that the fragmentation of quarks and gluons 
into K mesons is much weaker than the fragmentation 
into pions. Since h is the sum of charged pions and 
kaons, the inclusive charged pion cross section can be 
computed easily from this ratio. 

Similar conclusions can be drawn from the com- 
parisons shown in Fig. 6 and 7 with CDF data [5] at 
^/s = 630 GeV and y/s = 1.8 TeV, respectively. The 
definition of h is here again h = (h + + h~)/2. In 
both cases, the rapidity is averaged over the interval 
— 1 < y < 1. The data at both energies agree and are 
predicted best with scales equal to px- 

Our calculation seems to provide a rather good de- 
scription of data obtained in an energy range between 
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7. Same as in Fig. 6 for yfs = 1.8 TeV 
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and n — 10.64±0.31. We divide our theoretical predic- 
tions for the three choices of scales Pt/2, pt, and 2pT 
as well as the experimental data by the fit (4), taking 
for A, pt„ and n the central values. 




Pt (GeV) 



Fig. 8. The ratio theoretical cross section over fitted ansatz 
for the cross section as given by (4). The short-dashed, solid 
and long-dashed lines correspond to the NLO predictions 
from Fig. 2 for scales p = M = Mf set equal to Pt/2, pt 
and 2pT, compared to the ratio of experimental data points 
to fitted ansatz (4) from Fig. 2 (^/s = 500 GeV, \y\ < 2.5) 



200 GeV and 1.8 TeV, with an increase of energy by 
roughly one order of magnitude. Qualitatively our re- 
sults are independent of y/s at small pt, i.e. below 
3 GeV, and increase with increasing energy in the large 
Px-range, as it is characteristic for a hard-scattering 
cross section. 

The degree of agreement between data and theo- 
retical predictions, however, is rather difficult to assess 
by inspection of the logarithmic plots in Figs. 1-7, 
where the cross section drops rapidly over several or- 
ders of magnitude. It is clear that at small pr, i.e. 
Pt < 3 GeV, the theoretical curves deviate from the 
data. At higher pt, unfortunately, the data points have 
larger errors and establishing the level of agreement be- 
comes more problematic. 

To make our comparison easier, we employ a three 
parameter fit to the data of the form 



E 



fa 

d 3 p 



A(1+*L 

PT 



(4) 



In [3] the Vs = 500 GeV data of the UA1 Collabora- 
tion [3] have been fitted to (4) . This fit yields the pa- 
rameters A = 408±24mb/GcV 2 , p To = 1.61±0.08 GeV 
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Fig. 9. Same as Fig. 8 for the weaker gluon fragmentation 
given in [17] 

The result is shown in Fig. 8. It is clear that the 
fit reproduces the data rather well for pt < 3 GeV, 
but it is less satisfactory for pt around 3-4 GeV. It is 
compatible with the data for the largest values of Pt 
only thanks to the very large error of these last data 
points. The ratio of the central value for the data over 
the fit (4) seems to show an increase for increasing pt- 
This is the same behaviour shown by the theoretical 
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Fig. 10. Inclusive cross section for production of hadrons 
h = h + + hT as a function of pt for ^fs = 63 GeV and 
rapidity y = 0.75. Also shown are the data obtained by the 
CDHW Collaboration [2] at the ISR pp collider 

predictions for which the deviation from one becomes 
rather sharp at the largest values of pt here considered. 
It is more visible from this figure what was previously 
said, i.e. that, overall, the choice of scales <~ pt in our 
calculation seems to give the best agreement with the 
experimental data. 

It would be interesting to see in this same type 
of plots the effect of different fragmentation functions. 
We do not attempt this here, but we show in Fig. 9 
the same ratios obtained by using the parametrization 
of [17]. We remind here that this parametrization fits 
e + e~ data (except for the gluon fragmentation function 
from the OPAL measurements [16]) equally well as the 
parametrization of [9]. The gluon fragmentation func- 
tion in [17], however, is weaker than the one obtained 
in [9] where a fit to the OPAL data is also made. As 
one can see the results shown in Fig. 9 disagree with 
the data. The disagreement is stronger at smaller p^: 
the theoretical predictions are a factor 3-5 away from 
the data for pt < 3 GeV. A "strong" gluon fragmenta- 
tion to charged pions and kaons is therefore needed to 
explain the inclusive charged particle cross sections in 



So far all the inclusive particle production cross sec- 
tions from pp are for the sum of all charged particles 
only. No separation into particle species has been done. 
Such cross sections exist at lower center of mass ener- 
gies coming from the ISR pp collider [2] . These data are 
for charged pion production at y/s = 63 GeV, y ~ 0.75 
and a range of pr between 3 and 12 GeV. We have cal- 
culated this cross section with the same input as for 
the pp process (the p structure function in (2) is ob- 
viously replaced by the structure function for p) . The 
result is shown in Fig. 10 together with the data from 
the CDHW Collaboration [2] . The agreement with the 
theoretical curve with scales equal to px is excellent 
over the whole range of px- 

We have calculated inclusive single-charged hadron 
cross sections in full NLO, by using NLO structure 
functions, NLO fragmentation functions for charged 
pions and kaons and NLO hard scattering cross sec- 
tions. Our results were compared with experimental 
data from the CDHW, UAl, UA1-MIMI and CDF Col- 
laborations. We found very good agreement, in partic- 
ular with the UA1-MIMI data which have the smallest 
errors and extend over the largest pt range. The agree- 
ment with the data is satisfactory for all center of mass 
energies between 63 GeV and 1800 GeV in shape and 
absolute normalization, even in the small pt region. 
We have demonstrated that only those fragmentation 
functions with a large enough gluon contribution give a 
satisfactory account of the collider data. The strength 
of this gluon fragmentation function agrees with the 
OPAL data in the three-jet region sensitive to it [16]. 
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